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Abstract. The goal of these lectures, oriented towards the students just entering the field, is to 
provide an elementary introduction to QCD and the physics of nuclear interactions at high energies. 
We first introduce the general structure of QCD and discuss its main properties. Then we proceed to 
Glauber multiple scattering theory which lays the foundation for the theoretical treatment of nuclear 
interactions at high energies. We introduce the concept of Gribov's inelastic shadowing, crucial for 
the understanding of quantum formation effects. We outline the problems facing Glauber approach 
at high energies, and discuss how asymptotic freedom of QCD helps to resolve them, introducing 
the concepts of parton saturation and color glass condensate. 
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1. QUANTUM CHROMO-DYNAMICS - THE THEORY OF 
O : STRONG INTERACTIONS 

O 



1.1. What is QCD? 



Strong interaction is, indeed the strongest force of nature. It is responsible for over 
80% of the baryon masses, and thus for most of the mass of everything on Earth. Strong 
interactions bind nucleons in nuclei which, being then dressed with electrons and bound 
into molecules by the much weaker electro-magnetic force, give rise to the variety of the 
physical world. 

Quantum Chromodynamics (QCD) is the theory of strong interactions. The funda- 
mental degrees of freedom of QCD, quarks and gluons, are already well established 
even though they cannot be observed as free particles, but only in color neutral bound 
states (confinement). Today, QCD has firmly occupied its place as part of the Standard 
Model. However, understanding the physical world does not only mean understanding its 
fundamental constituents; it means mostly understanding how these constituents interact 
and bring into existence the entire variety of physical objects composing the universe. In 
these lectures, we try to explain why high energy nuclear physics offers us unique tools 
to study QCD. 



7.7.7. The QCD Lagrangian 

So what is QCD? QCD emerges when the naive quark model is combined with local 
SU(3) gauge invariance. Quark model classifies the large number of hadrons in terms 



of a few, more fundamental constituents. Baryons consist of three quarks, while mesons 
are made of a quark and an antiquark. For example, the proton is made of two up-quarks 
and one down quark, \p) = \uud), and the 7T + -meson contains one up and one anti- 
down quark, \ k + ) = \ud). However, the quark model in this naive form is not complete, 
because the Pauli exclusion principle would not allow for a particle like the A isobar 
|A ++ ) = \uuu) with spin 3/2. The only way to construct a completely antisymmetric 
wavefunction for the A ++ is to postulate an additional quantum number, which may 
be called "color". Quarks can then exist in three different color states; one may choose 
calling them red, green and blue. Correspondingly, we can define a quark-state "vector" 
with three components, 

/ q Kd (x) \ 

q{x) = <? green W • (1) 

V <? blue to / 

The transition from quark model to QCD is made when one decides to treat color simi- 
larly to the electric charge in electrodynamics. As is well known, the entire structure of 
electrodynamics emerges from the requirement of local gauge invariance, i.e. invariance 
with respect to the phase rotation of electron field, exp(ia(x)), where the phase a de- 
pends on the space-time coordinate. One can demand similar invariance for the quark 
fields, keeping in mind that while there is only one electric charge in QED, there are 
three color charges in QCD. 

To implement this program, let us require the free quark Lagrangian, 

^free= E E 9 to J~ ~ ™q ) 1 to , ( 2 ) 

q=u,d,s... colors \ H / 

to be invariant under rotations of the quark fields in color space, 

U: q j (x) -> U jk (x)q k (x), (3) 

with j,k E {1 . . .3} (we always sum over repeated indices). Since the theory we build 
in this way is invariant with respect to these "gauge" transformations, all physically 
meaningful quantities must be gauge invariant. 

In electrodynamics, there is only one electric charge, and gauge transformation in- 
volves a single phase factor, U = exp{ia(x)). In QCD, we have three different colors, 
and U becomes a (complex valued) unitary 3x3 matrix, i.e. U^U = 1/1/^ = 1, with de- 
terminant Det U = 1 . These matrices form the fundamental representation of the group 
SU (3) where 3 is the number of colors, N c = 3. The matrix U has N% — 1 = 8 inde- 
pendent elements and can therefore be parameterized in terms of the 8 generators T^-, 
a G { 1 ... 8} of the fundamental representation of SU(3), 

C/(x) = exp(-i0 fl (x)r) (4) 

By considering a transformation U that is infinitesimally close to the 1 element of the 
group, it is easy to see that the matrices T a must be Hermitian (T a = T a ^) and traceless 
(tr T a = 0). The T a, s do not commute; instead one defines the SU (3) structure constants 



fabc tne commutator 

' abc 



if» h J c . (5) 



These commutator terms have no analog in QED which is based on the abelian gauge 
group 2/(1). QCD is based on a non-abelian gauge group SU(3) and is thus called a 
non-abelian gauge theory. 

The generators T a are normalized to 

ixT a T b = l -8 ab , (6) 

where 8 ab is the Kronecker symbol. Useful information about the algebra of color 
matrices, and their explicit representations, can be found in many textbooks (see, e.g., 

Since U is ^-dependent, the free quark Lagrangian (fj) is not invariant under the 
transformation @. In order to preserve gauge invariance, one has to introduce, following 
the familiar case of electrodynamics, the gauge (or "gluon") field A^. (x) and replace the 

derivative in @ with the so-called covariant derivative, 

aVto - D^(x)^{8 kj d^-iA^x)} q \x). (7) 

Note that the gauge field A£. (x) = A^T?-(x) as well as the covariant derivative are 3x3 

matrices in color space. Note also that Eq. Q7|) differs from the definition often given 
in textbooks, because we have absorbed the strong coupling constant in the field AP. 
With the replacement given by Eq. (0), all changes to the Lagrangian under gauge 
transformations cancel, provided A^ transforms as 

U : A^(jc) -> U(x)A^(x)U\x) +iU(x)d^U\x). (8) 

(From now on, we will often not write the color indices explicitly.) 
The QCD Lagrangian then reads 

^ QCD = Y,q(x){iY^-m q )q(x)-^- I trG^(x)G^(x) , (9) 



vQQQgaaaasa 




FIGURE 1. Due to the non-abelian nature of QCD, gluons carry color charge and can therefore interact 
with each others via these vertices. 



where the first term describes the dynamics of quarks and their couplings to gluons, 
while the second term describes the dynamics of the gluon field. The strong coupling 
constant g is the QCD analog of the elementary electric charge e in QED. The gluon 
field strength tensor is given by 

G» v (x) =i[D^,D v ] = d^A v (x)-d v A^(x)~i[A^(x),A v (x)]. (10) 

This can also be written in terms of the color components A^ of the gauge field, 

GTM = d^Alix) -d v A£( X ) +f abc Al{x)A v c {x). (11) 

For a more complete presentation, see [0] and modern textbooks like [Jl|, |3|, 

The crucial, as will become clear soon, difference between electrodynamics and QCD 
is the presence of the commutator on the r.h.s. of Eq. (flOl). This commutator gives rise 
to the gluon-gluon interactions shown in Fig. [1] that make the QCD field equations non- 
linear: the color fields do not simply add like in electrodynamics. These non-linearities 
give rise to rich and non-trivial dynamics of strong interactions. 



7.7.2. Asymptotic Freedom 

Let us now turn to the discussion of the dynamical properties of QCD. To understand 
the dynamics of a field theory, one necessarily has to understand how the coupling 
constant behaves as a function of distance. This behavior, in turn, is determined by the 
response of the vacuum to the presence of external charge. The vacuum is the ground 
state of the theory; however, quantum mechanics tells us that the "vacuum" is far from 
being empty - the uncertainty principle allows particle-antiparticle pairs to be present in 
the vacuum for a period time inversely proportional to their energy. In QED, the electron- 
positron pairs have the effect of screening the electric charge, see Fig. 0. Thus, the 
electromagnetic coupling constant increases toward shorter distances. The dependence 
of the charge on distance is given by 

eHr) = ffi' , (12) 

which can be obtained by resumming (logarithmically divergent, and regularized at the 
distance r ) electron-positron loops dressing the virtual photon propagator. 

The formula ( |12| ) has two surprising properties: first, at large distances r away from the 
charge which is localized at r Q , r ^> r Q , where one can neglect unity in the denominator, 
the "dressed" charge e(r) becomes independent of the value of the "bare" charge e(r Q ) - 
it does not matter what the value of the charge at short distances is. Second, in the local 
limit r Q — > 0, if we require the bare charge e(r ) be finite, the effective charge vanishes 
at any finite distance away from the bare charge! This is the celebrated Landau's zero 
charge problem [f5p: the screening of the charge in QED does not allow to reconcile 
the presence of interactions with the local limit of the theory. This is a fundamental 
problem of QED, which shows that i) either it is not a truly fundamental theory, or ii) 



V 

FIGURE 2. /n QED, virtual electron-positron pairs from the vacuum screen the bare charge of the elec- 
tron. The larger the distance, the more pairs are present to screen the bare charge and the electromagnetic 
coupling decreases. Conversely, the coupling is larger when probed at short distances. 

Eq. (|T2|), based on perturbation theory, in the strong coupling regime gets replaced by 
some other expression with a more acceptable behavior. The latter possibility is quite 
likely since at short distances the electric charge becomes very large and its interactions 
with electron-positron vacuum cannot be treated perturbatively. A solution of the zero 
charge problem, based on considering the rearrangement of the vacuum in the presence 
of "super-critical", at short distances, charge was suggested by Gribov [§]. 

Fortunately, because of the smallness of the physical coupling a em (r) = e 2 (r) / (4tt) = 
1/137, this fundamental problem of the theory manifests itself only at very short dis- 
tances ~ exp(— 3/[8a em ]). Such short distances will probably always remain beyond the 
reach of experiment, and one can safely apply QED as a truly effective theory. 

In QCD, as we are now going to discuss, the situation is qualitatively different, and 
corresponds to anrf- screening - the charge is small at short distances and grows at larger 
distances. This property of the theory, discovered by Gross, Wilczek, and Politzer [[7]], is 
called asymptotic freedom. 

While the derivation of the running coupling is conventionally performed by using 
field theoretical perturbation theory, it is instructive to see how these results can be 
illustrated by using the methods of condensed matter physics. Indeed, let us consider 
the vacuum as a continuous medium with a dielectric constant £. The dielectric constant 
is linked to the magnetic permeability /i and the speed of light c by the relation 

£H = \ = l. (13) 

cr 

Thus, a screening medium (e > 1) will be diamagnetic (/x < 1), and conversely a 
paramagnetic medium (fi > 1) will exhibit antiscreening which leads to asymptotic 
freedom. In order to calculate the running coupling constant, one has to calculate the 
magnetic permeability of the vacuum. We follow [|]] in our discussion, where this has 



been done in a framework very similar to Landau's theory of the diamagnetic properties 
of a free electron gas. In QED one has 



2e (r ) r 

%™=l + ^ln->l 04) 

So why is the QCD vacuum paramagnetic while the QED vacuum is diamagnetic? The 
energy density of a medium in the presence of an external magnetic field B is given by 

u = - 1 -AtzxB 2 (15) 

where the magnetic susceptibility % is defined by the relation 

\l = \+Ail%. (16) 

When electrons move in an external magnetic field, two competing effects determine the 
sign of magnetic susceptibility: 

• The electrons in magnetic field move along quantized orbits, so-called Landau 
levels. The current originating from this movement produces a magnetic field with 
opposite direction to the external field. This is the diamagnetic response, % < 0. 

• The electron spins align along the direction of the external 5-field, leading to a 
paramagnetic response (x > 0). 

In QED, the diamagnetic effect is stronger, so the vacuum is screening the bare charges. 
In QCD, however, gluons carry color charge. Since they have a larger spin (spin 1) 
than quarks (or electrons), the paramagnetic effect dominates and the vacuum is anti- 
screening. 

Let us explain this in more detail. Basing on the considerations given above, the 
energy density of the QCD vacuum in the presence of an external color-magnetic field 
can be calculated by using the standard formulas of quantum mechanics, see e.g. 0], 
by summing over Landau levels and taking account of the fact that gluons and quarks 
give contributions of different sign. Note that a summation over all Landau levels would 
lead to an infinite result for the energy density. In order to avoid this divergence, one 
has to introduce a cutoff A with dimension of mass. Only field modes with wavelength 
X > 1/A are taken into account. The upper limit for A is given by the radius of the 
largest Landau orbit, r Q ~ 1 / y/gB, which is the only dimensionful scale in the problem; 
the summation thus is made over the wave lengths satisfying 

1 „ 1 

^A>-, (17) 



The result is [§] 

nrn i AlN c -2N f , a 2 

M — -2* 2 ^8^ 2ln W (18) 

where Nj> is the number of quark flavors, and N c = 3 is the number of flavors. Comparing 
this with Eqs. ([15]) and (|16|), one can read off the magnetic permeability of the QCD 



vacuum, 

UN c -2N f , a 2 

«) = i + ^r^rm M >i. (19) 

The first term in the denominator (11/V C ) is the gluon contribution to the magnetic 
permeability. This term dominates over the quark contribution (2A^) as long as the 
number of flavors N* is less than 17 and is responsible for asymptotic freedom. 
The dielectric constant as a function of distance r is then given by 



(20) 

\gB\->l/r 



The replacement \J\gB\ — > 1/r follows from the fact that e and /i in Eq. ( pH ) should be 
calculated from the same field modes: the dielectric constant e(r) could be calculated by 
computing the vacuum energy in the presence of two static colored test particles located 
at a distance r from each other. In this case, the maximum wavelength of field modes 
that can contribute is of order r so that 

r^A^i (21) 
Combining Eqs. ([17|) and (|2l)), we identify r = 1/ \J\gB\ and find 

= nNc -il f ~ - - < L ( 22 ) 




FIGURE 3. r/;e running coupling constant OC s (Q 2 ) as a function of momentum transfer Q 2 determined 
from a variety of processes. The figure is from /|7C|/, courtesy ofS. Bethke. 



With a s (r l )/a s (r 2 ) = £@£ D (r 2 ) /££c D ( r i) one finds t0 lowest order in a s 

°^l)= UN c -2N f 2 : /.y ( 23 ) 

Apparently, if r\ < r 2 then & s {r\) < a s(r 2 ). The running of the coupling constant is 
shown in Fig. ^ Q ~ 1/r. The intuitive derivation given above illustrates the original 
field-theoretical result of []7p. 

At high momentum transfer, corresponding to short distances, the coupling constant 
thus becomes small and one can apply perturbation theory, see Fig. [3|. There is a va- 
riety of processes that involve high momentum scales, e.g. deep inelastic scattering, 
Drell-Yan dilepton production, e + e -annihilation into hadrons, production of heavy 
quarks/quarkonia, high p T hadron production QCD correctly predicts the Q 2 de- 
pendence of these, so-called "hard" processes, which is a great success of the theory. 



1.2. Challenges in QCD 

7.2.7. Confinement 

While asymptotic freedom implies that the theory becomes simple and treatable at 
short distances, it also tells us that at large distances the coupling becomes very strong. 
In this regime we have no reason to believe in perturbation theory. In QED, as we have 
discussed above, the strong coupling regime starts at extremely short distances beyond 
the reach of current experiments - and this makes the "zero-charge" problem somewhat 
academic. In QCD, the entire physical World around us is defined by the properties of 
the theory in the strong coupling regime - and we have to construct accelerators to study 
it in the much more simple, "QED-like", weak coupling limit. 

We do not have to look far to find the striking differences between the properties 
of QCD at short and large distances: the elementary building blocks of QCD - the 
"fundamental" fields appearing in the Lagrangean @, quarks and gluons, do not exist 
in the physical spectrum as asymptotic states. For some, still unknown to us, reason, 
all physical states with finite energy appear to be color-singlet combinations of quarks 
and gluons, which are thus always "confined" at rather short distances on the order of 1 
fm. This prevents us, at least in principle, from using well-developed formal S-matrix 
approaches based on analyticity and unitarity to describe quark and gluon interactions. 

The property of confinement can be explored by looking at the propagation of heavy 
quark-antiquark pair at a distance R propagating in time a distance T. An object which 
describes the behavior of this system is the Wilson loop [ pT] ] 

W(R,T) = Tr 

where A a ^ is the gluon field, T a is the generator of SU (3), and the contour C is chosen 
as a rectangle with side R in one of the space dimensions and T in the time direction. It 



P exp 



i J A^T'dx? 



(24) 



can be shown that at large T the asymptotics of the Wilson loop is 

lim W(R, T) = exp [—TV(R)] , (25) 

T — >oo 

where V(R) is the static potential acting between the heavy quarks. At large distances, 
this potential grows as 

V(R) = oR, (26) 

where o ~ 1 GeV/fm is the string tension. We thus conclude that at large T and R the 
Wilson loop should behave as 

W(R,T) ~exp[-oTR], (27) 



The formula ( |27| ) is the celebrated "area law", which signals confinement. 

It should be noted, however, that the introduction of dynamical quarks leads to the 
string break-up at large distances, and the potential V(R) saturates at a constant. The 
presence of light dynamical quarks is most important in Gribov's confinement scenario 
|f|], in which the color charges at large distances behave similarly to the "supercriti- 
cal" charge in electrodynamics, polarizing the vacuum and producing copious quark- 
antiquark pairs which screen them. In this scenario, in the physical world with light 
quarks there is never a confining force acting on color charges at large distances, just 
quark-antiquark pair production ("soft confinement"). This may explain why the spec- 
tra of jets, for example, computed in perturbative QCD, appear to be consistent with 
experiment; this fact would be difficult to reconcile with the existence of strong con- 



fining forces. There exists a special situation, however, when the law fl27| ) should be 
appropriate even in the presence of light quarks - the heavy quarkonium. The sizes of 
heavy quarkonia are quite small, and their masses are below the threshold to produce 
a pair of heavy mesons. This is why heavy quarkonia are especially useful probes of 
confinement. 

At high temperatures, the long-range interactions responsible for confinement be- 
come screened away - instead of the growing potential (^), we expect 

V(R) ~ -^p- exp(-m D R), (28) 
K 

where m D ~ gT is the Debye mass. Mathematically, this transition to the deconfined 
phase can again be studied by looking at the properties of the Wilson loop. At finite 
temperature, the theory is defined on a cylinder: Euclidean time % varies within < T < 
/3 = 1/r, and the gluon fields satisfy the periodic boundary conditions: 

A;&0)=A£(*,jB). (29) 

Let us now consider the Wilson loop wrapped around this cylinder (the Polyakov loop), 
and choose a gauge where Aq is time-independent: 

P(x) = Tr exp [igjBAg(Jc )t a ] ; (30) 

the correlation function of these objects can be defined as 



C T (x) =<P(x)P*(x) > T 



(31) 



Again, it can be shown that this correlation function is related to the free energy, and 
thus static potential V(R), of the heavy quark-antiquark pair. Assuming, as before, that 
the heavy quarks are separated by the spatial distance R = \x\, one finds 

C T (R)~exp[-pV(R)]. (32) 

Again, if we define the limit value L{T) of the correlation function, 

lim C T (R)=L(T) (33) 

R— >oo 

it would have to vanish in the confined phase in the absence of dynamical quarks, since 
V(R) tends to infinity in this case: L(T) = 0. In the deconfined phase, on the other hand, 
because of the screening V(R) should tend to a constant, and this implies a finite value 
L(T) ^ 0. The correlation function of Polyakov loops therefore can be used as an order 
parameter of the deconfinement. The behavior of L(T) as a function of temperature has 
been measured on the lattice; one indeed observes a transition from the confined phase 
with L(T) = to the deconfined phase with L(T) ^ at some critical temperature T c . 
In the presence of light quarks, as we have already discussed above, the potential would 
tend to a constant even in the confined phase, and L(T) ceases to be a rigorous order 
parameter. 



7.2.2. Chiral symmetry breaking 

The decades of experience with "soft pion" techniques and current algebra convinced 
physicists that the properties of the world with massless pions are quite close to the 
properties of our physical World. The existence of massless particles is always a mani- 
festation of a symmetry of the theory - photons, for example, appear as a consequence 
of local gauge invariance of the electrodynamics. However, unlike photons, pions have 
zero spin and cannot be gauge bosons of any symmetry. The other possibility is pro- 
vided by the Goldstone theorem, which states that the appearance of massless modes 
in the spectrum can also reflect a spontaneously broken symmetry, i.e. the symmetry of 
the theory which is broken in the ground state. Because of the great importance of this 
theorem, let us briefly sketch its proof. 

Suppose that the Hamiltonian H of the theory is invariant under some symmetry 
generated by operators Q t , so that 

[H,Q f ]=0. (34) 

Spontaneous symmetry breaking in the ground state of theory implies that for some of 
the generators Q i 

G/|0>^0. (35) 

Since Q i commute with the Hamiltonian, this means that this new state <2-|0 > has the 
same energy as the ground state. The vacuum is therefore degenerate, and in a rela- 
tivistically invariant theory this implies the existence of massless particles - Goldstone 



bosons. A useful example of that is provided by the phonons in a crystal, where the con- 
tinuous translational symmetry of the QED Lagrangean is spontaneously broken by the 
existence of the fixed period of the crystal lattice. 

Even though all six quark flavors enter the Lagrangean, it is intuitively clear that at 
small scales Q << M c ,M,,M t , heavy quarks should not have any influence on the dy- 
namics. In a rigorous way this statement is formulated in terms of decoupling theorems, 
which we will discuss in detail later. At the moment let us just assume that we are inter- 
ested in the low-energy behavior, and that only light quarks are relevant for that purpose. 
Then it makes sense to consider the approximate symmetry, which becomes exact when 
the quarks are massless. In fact, in this limit, the Lagrangean does not contain any terms 
which connect the right- and left-handed components of the quark fields: 

1 1 

<lR=2^ + Y 5 )q\ ^ L = 2^~ r ^ q - (36) 

The Lagrangean of QCD @ is therefore invariant under the independent transformations 
of right- and left-handed fields ("chiral rotations"). In the limit of massless quarks, QCD 
thus possesses an additional symmetry U L {Nj) x U R (N^) with respect to the independent 

transformation of left- and right-handed quark fields q L R = j( 1 ± y 5 )g: 

9 L ^V L q L ; 1r^Vr<1r> V L ,V R eU(N f ); (37) 

this means that left- and right-handed quarks are not correlated. 

Even a brief look into the Particle Data tables, or simply in the mirror, can convince 
anyone that there is no symmetry between left and right in the physical World. One thus 
has to assume that the symmetry ( |3"7| ) is spontaneously broken in the vacuum. 

The presence of the "quark condensate" < qq > in QCD vacuum signals spontaneous 
breakdown of this symmetry, since 

<qq>=<q~ L q R > + <q~ R q L >, (38) 

which means that left- and right-handed quarks and antiquarks can transform into 
each other. Quark condensate therefore can be used as an order parameter of chiral 
symmetry. Lattice calculations show that around the deconfinement phase transition, 
quark condensate dramatically decreases, signaling the onset of the chiral symmetry 
restoration. 

This spontaneous breaking of U L (3) x U R (3) chiral symmetry, by virtue of the Gold- 
stone theorem presented above, should give rise to 3 2 = 9 Goldstone particles. The flavor 
composition of the existing eight candidates for this role (3 pions, 4 kaons, and the 77) 
suggests that the U A (1) part of U L {3) xU R (3) =SU L (3) xSU R (3) xU v {l) *U A (1) does 
not exist. This constitutes the famous "U A (1) problem". 

1.2.3. The origin of mass 

There is yet another problem with the chiral limit in QCD. Indeed, as the quark masses 
are put to zero, the Lagrangian (Q) does not contain a single dimensionful scale - the 




FIGURE 4. The triangle graph that leads to the U A {\)- anomaly. The corresponding graph with the two 
gluons interchanged in the final state is not shown. 



only parameters are pure numbers N c and Nj. The theory is thus apparently invariant 
with respect to scale transformations, and the corresponding scale current is conserved: 
dfi s fi = 0- However, the absence of a mass scale would imply that all physical states in 
the theory should be massless! 



1.2.4. Quantum anomalies 

Both apparent problems - the missing U A (l) symmetry and the origin of hadron 
masses - are related to quantum anomalies. A symmetry of a classical theory can be 
broken when that theory is quantized, due to the requirements of regularization and 
renormalization. This is called anomalous symmetry breaking. Regularization of the 
theory on the quantum level brings in a dimensionful parameter - remember the cutoff 
A of Eq. ( |T7| ) we had to impose on the wavelength of quarks and gluons. 

Once the theory is quantized, we already know that the coupling constant is scale 
dependent and therefore scale invariance is broken (note that the four-divergence of the 
scale current in field theory is equal to the trace of the energy momentum tensor ©f,). 
One finds 

d% = % = 2> 9 «?<7+ ^W V G MV , (39) 
where /3 (g) is the QCD /3-function, which governs the behavior of the running coupling: 

^ = /■(,): 

note that as discussed in Section 1.1.1 we include coupling g in the definition of the 
gluon fields. As we already discussed, at small coupling g, the /3 function is negative, 
which means that the theory is asymptotically free. The leading term in the perturbative 
expansion is (compare with Eq. (E31)) 



P(g) = -b-^ + 0(g 5 ), b = UN c -2N f , (41) 

where N c and Nj- are the numbers of colors and flavors, respectively. 

Hadron masses are related to the forward matrix element of trace of the QCD energy- 
momentum tensor, 2mj t = (h\&^ \h). Apparently, light hadron masses must receive 



dominant contributions from the G 2 -term in Eq. (f39|). Note also that the flavor sum 
in Eq. (T33|) includes heavy flavors, too. This would lead to the unphysical picture 
that e.g. the proton mass is dominated by heavy quark masses. However, the heavy 
flavor contribution to the sum is exactly canceled by a corresponding heavy flavor 



contribution to the /3 -function. 

Similar thing happens with the axial current, & = qy^y 5 q, generated by the U A (1) 
group. The corresponding axial charge is not conserved because of the contribution of 
the triangle graph in Fig. |], and the four-divergence of the axial current is given by [ |T2"| ] 



= Y,2im q qr 5 q + ^tTG^G^ v , (42) 

where G^ v = e u VK iG K ^ /2 is the dual field strength tensor. Since the gluonic part on 
the rhs of this equation is a surface term (a full divergence), there would be no physical 
effect, if the QCD vacuum were "empty". 



1.2.5. Classical solutions 

However, it appears that due to non-trivial topology of the SU (3) gauge group, 
QCD equations of motion allow classical solutions even in the absence of external 
color source, i.e. in the vacuum. The well-known example of a classical solution is 
the instanton, corresponding to the mapping of a three-dimensional sphere S 3 onto the 
SU (2) subgroup of color SU (3) (for reviews, see [ |13| , [14|]). As a result, the ground state 
of classical Chromodynamics is not unique. There is an enumerable infinite number of 
gauge field configurations with different topologies (corresponding to different winding 
number in the S 3 — > SU (2) mapping), and the ground state looks like a periodic potential, 
see Fig. |5| 

In a quantum theory, however, the system will not stay in one of the minima, like 
the classical system would. Instead, there will be tunneling processes between different 
minima. These tunneling processes, in Minkowski space, correspond to instantons. Since 
tunneling, in general, lowers the ground state energy of the system, one expects the QCD 
vacuum to have a complicated structure. 

Instantons, through the anomaly relation (R2I), lead to the explicit violation of the 



U A (l) symmetry and thus solve the mystery of the missing ninth Goldstone boson - the 
r\' . Physically, axial symmetry UA\) is broken because the tunneling processes between 
topologically different vacua are accompanied by the change in quark helicity - even in 
the vacuum, left-handed quarks periodically turn into right-handed and vice versa. 



1.2.6. Strong CP problem 

The vacuum structure shown in Fig. |5] immediately leads to a puzzle known as the 
strong CP problem: When one calculates the expectation value of an observable in the 
vacuum, one has to average over all topological sectors of the vacuum. This is equivalent 



A A A A 

^(0) ^(1) ^(2) ^(3) 

FIGURE 5. Topological structure of QCD vacuum. The minima correspond to classical ground states 
with topologically different gauge field configurations A^ t y Also shown is an instanton trajectory interpo- 
lating between the classical vacua A,^. and A^. The third axis shows the Euclidean time x. From /[T^j/; 
courtesy ofH. Forkel. 



to adding an additional term to the QCD-Lagrangian, 

~ 

■^QCD ~^ =^QCD — 16^2 tr G^G^y, (43) 

where G [0, lit] is a parameter of the theory which has to be determined from exper- 
iment. Since the 0-term in Eq. is CP violating, a non-zero value of would have 
immediate phenomenological consequences, e.g. an electric dipole moment of the neu- 
tron. However, precision measurements of this dipole moment constrain to < 10 -9 . 
The fact that is so unnaturally small constitutes the strong CP problem. The most 
likely solution to this problem [ ]T3| ] implies the existence of a light pseudoscalar meson, 
the axion. However, despite many efforts, axions remain unobserved in experiment. 



7.2.7. Phase structure 

As was repeatedly stated above, the most important problem facing us in the study 
of all aspects of QCD is understanding the structure of the vacuum, which, in a manner 
of saying, does not at all behave as an empty space, but as a physical entity with a 
complicated structure. As such, the vacuum can be excited, altered and modified in 
physical processes 

Collisions of heavy ions are the best way to create high energy density in a "macro- 
scopic" (on the scale of a single hadron) volume. It thus could be possible to create and 



to study a new state of matter, the Quark-Gluon Plasma{QGP), in which quarks and 
gluons are no longer confined in hadrons, but can propagate freely. The search for QGP 
is one of the main motivations for the heavy ion research. 

Lattice calculations predict that QCD at high temperatures undergoes phase transi- 
tions in which confinement property is lost and chiral symmetry is restored. The critical 
temperature for the chiral phase transition is similar (or maybe even equal) to the critical 
temperature for deconfinement. 

Heavy ion collisions at RHIC may also give us the possibility to study the 6 angle 
dependence of the QCD phase diagram. In a heavy ion collision, bubbles containing a 
metastable vacuum with 0^0 may be produced, and reveal themselves through their 
unusual decay pattern []T7|]. 

2. NUCLEAR INTERACTIONS AT HIGH ENERGIES 

2.1. Glauber-Gribov Theory 

It is intuitively clear that heavy ion collisions are governed by multiple scattering 
effects. As a short introduction to the basics of multiple scattering theory, we introduce 
here the eikonal approximation to high energy scattering processes and the Glauber 
multiple scattering theory [JTSp. We also discuss Gribov's inelastic corrections [ |TP| ] to 
Glauber's theory. 



2.7.7. The Eikonal Approximation 

The eikonal approximation is the classical approximation to the angular momentum 
/. In partial wave expansion, i.e. in an expansion in angular momentum eigenstates, the 
scattering amplitude f(s,t) reads [j^] 

P z (cos0), (44) 

where s and t are the usual Mandelstam variables (center-of-mass energy squared and 
invariant momentum transfer, respectively), p is the momentum of the projectile and P, 
are the Legendre functions, which depend on the cosine of the scattering angle 9. All 
information about the interaction is contained in the scattering phases 8 r 

High energy scattering is of course a process that is far from being spherically 
symmetric. Therefore, very large values of / will dominate the sum Eq. ( f44l ) and we 
can treat the angular momentum classically. Since the angular momentum is given by 
pb, one replaces the variable / by the impact parameter b, 

pb = l+ 1 -. (45) 
Note that b is now a continuous variable, so angular momentum is no longer quantized. 



fM = — 1(21 + 1) 



e 2i5, 



At large I and for small scattering angles 6, the Legendre functions can be expressed 
to good approximation as 
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(46) 



where q = p — p' is the momentum transfer in the scattering process (t « — \q\) and 
|p| = for elastic scattering. At high energy, q lies in the impact parameter plane. We 
have used the relation 



(2Z + l)sin(0/2)cos(0) = 2psin(0/2)^-^cos(0) =qb 

P 

to obtain the second equality in Eq. (|46|). 

Thus, the scattering amplitude in eikonal approximation reads 



(47) 
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where the phase shift of the projectile is related to the scattering phase 8 l by 

X(s,b)=25(s,b). 

In the case of scattering off a potential V(r), this phase shift is simply given by 

1 



Zip) 



V(f)dz, 



where v is the velocity of the projectile. The scattering amplitude then reads 



f(s,t) = ^Jd 2 bj* b 



1 — exp 



V{r)dz 



(48) 



(49) 



(50) 



(51) 



The total cross section can now be obtained from the forward scattering amplitude via 
the optical theorem, 



a tot = —Jmf{s,t = 0)=2j d 2 b (l -Ree^W) . 



(52) 



For completeness, we also give the expressions for the elastic and inelastic cross 
sections. The elastic cross section is obtained by squaring the elastic scattering amplitude 
and integrating over the solid angle, 



a el = j d£l pl \f(dM 2 - 



(53) 
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predominantly in forward direction, one obtains 



With the approximation dQ. p , ~ d 2 p' /p' 2 , which assumes that scattering takes place 
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FIGURE 6. Scattering off a composite system. The impact parameter of the projectile is denoted by b, 
while the impact parameters of the scattering centers are denoted by Sj. 



Finally, the inelastic cross section is 



a inel = °tot - Orf = / d 2 b ( 1 



Mb) 



(55) 



For potential scattering, the inelastic cross section, of course, vanishes because %{b) is 
real. In general, however, %(b) will have an imaginary part. 

The expressions Eqs. d52|), ( |54l ) and ( |55|) could have been obtained directly from the 
partial wave decomposition of the total, elastic and inelastic cross section, as well. The 
conditions under which the eikonal approximation is applicable are investigated in detail 
in [ph. 



2.7.2. Multiple Scattering Theory 

Based on the eikonal approximation, it is quite straightforward to develop a theory 
for scattering off a composite system. In this section, we explain the basic features 
of the multiple scattering theory developed by Glauber [|l~8]] ■ A much more detailed 
presentation of this subject can be found in [jT8p. 

Assuming that the scatterings on different nucleons are independent, the phase shifts 
from each scattering simply add up, 

A 

X(b, r v r 2 ,... r a ) = £ Xj@ ~ ?y) • (56) 

;=t 

Here, b is the impact parameter of the projectile and j = I... A are the impact 
parameters of the A nucleons in the nucleus, see Fig. The amplitude for scattering 



off a nuclear target then can be written as 



Ffi = ^jd%^ r \f\\-f[^-^\f) (57) 



2n 



d^^ifii-Yiii-rjib-sjm, (58) 



where |/) and |i) are the final and initial state of the target, respectively. In the second 
step, we introduced the profile function y(b), which is related to the single-scattering 
amplitude f(q) by 

7{b) = ^-Jd 2 q^ rb f{q). (59) 

Thus, we have expressed the nuclear scattering amplitude in terms of the amplitude for 
scattering off a single nucleon. 

In the case of a purely imaginary f(q), y(b) is the probability of absorption of the 
projectile by a nucleon and the nuclear scattering amplitude, Eq. ( |58| ) has a simple prob- 
abilistic interpretation. Namely, 1 — y-(b — Sj) is the probability of not being absorbed 
by nucleon number j. Taking the product over all j 6 {1 .. .A] yields the probability of 
not being absorbed by any nucleon in the target. Finally, 1 — YYj=\ [1 ~~ 7j(b ~~ */)] is me 
probability that the projectile is absorbed by any of the nucleons. 

Also, if one in addition assumes that all nucleons in the target are identical, the nuclear 
cross section can be expressed in terms of the cross section for scattering on a single 
nucleon, 

att = — lmftf(f = 0) (60) 



P 



^( l .( 1 .-jm ] , (61) 

d 2 bl l-^J- a "» T * (b) I I ,62. 



where the nuclear thickness function TAb) is the integral over the nuclear density, 

T A {b)= / dzp A (b,z). (63) 




The simple expression, Eq. (|6T|), resums all multiple scattering terms. We stress that 
the probabilistic interpretation of Eq. ( |5B| ) as well as Eqs. (|^T|) and (]62|) only hold for a 
purely imaginary f(q) . 

The meaning of the nuclear scattering amplitude, Eq. (|5%|), is further explained by 
expanding the probability of particle absorption by any of the nucleons in powers of 



FIGURE 7. Illustration of the single and double scattering terms in Eq. (p3[). The coherent sum over 
all graphs leads to interferences that reduce the total cross section. 



A 



(64) 



7=1 



A 



E Y k (b ~ s k ) - £ Y k (b - ? k )y.(S -?.) + ■■■■ (65) 



The first two terms in this expansion are illustrated in Fig. [7|. The first term in Eq. ( |65| ) is 
just the sum of single scattering amplitudes. However, different nucleons in the nucleus 
compete to interact with the projectile. This effect is contained in the second term in 
Eq. (|63|), which reduces the cross section. This reduction is an interference effect that 
appears because the amplitudes for scattering on different nucleons have to be added 
coherently. This destructive interference can be observed in experiment as shadowing in 
hadron-nucleus interactions (eclipse effect in deuterium). Note, however, that shadowing 
is not completely explained by Glauber theory, as will be explained in the following 
section. 

The easiest application of Glauber multiple scattering theory to nuclear systems is the 
calculation of the inelastic nucleus-nucleus (AB) cross section, which can be written as 




(66) 



here, P Q (b) is the probability that no interaction takes place, 



W = d-«)) 



.AB 



(67) 



where the nuclear overlap function is given by 




(68) 



(Obviously, 1 —P Q (b) is then the probability of an inelastic interaction, and the meaning 
of Eq. ( |66| ) becomes very transparent.) As it is common, we have labeled the two nuclei 
by their atomic mass numbers A and B. 



Another application is the calculation of inclusive particle spectra. With the help 
of crossing symmetry, the cross section for production of a particle of type a in an 
AB collision, AB — > aX, can be calculated from the total cross section of the process 
dAB — > X, where a is the antiparticle of a. According to so-called AGK cutting rules 
[|2l|], the nuclear cross section for this process is given by 

d 3 n a d 3 a a 
17 AB -Tab&E-^- (69) 



d 2 bd 3 p ASW d 3 p 
Integration over impact parameter b yields 



E ^k =AB E^L, (70) 
d i p d i p 

and correspondingly the charged particle multiplicity would scale proportional to AB, 

d^ 1=AB J_d_^ k 

dJ] dJ] 

the meaning of which is obvious - if collisions are truly independent, the resulting 
multiplicity should scale with the number of collisions, AB. 

However, the relation ( [77] ) appears to be badly violated in experiment. What went 
wrong? It appears that the disagreement between the result Eq. (jTTJ) and experimental 
data is due to the fact that there are important corrections to the Glauber multiple 
scattering theory, which we neglected so far. These corrections are known as Gribov's 
inelastic shadowing []T9]] and will be the subject of the next section. 



2.1.3. Gribov's "Inelastic Shadowing" 

The assumed independence of nucleon-nucleon collisions is violated by the diagrams 
of the type of Fig. |[ where the projectile is excited into a state \n) by the interaction. 
The diagram in Fig. |]does not describe independent collisions, and at high energies it 
will interfere with the double scattering graph in Fig. [7J. 

The excitation of an inelastic state in the scattering is accompanied by a longitudinal 
momentum transfer 

where Mr is the invariant mass of the excited system and M i is the invariant mass of 
the projectile in the initial state. The diagram in Fig. [8| is only important if it can make 
a significant contribution to the forward scattering amplitude F». This requires that the 
longitudinal momentum transfer must be so small that the nucleus has a chance to remain 
intact, i.e. 

*P L R A &1, (73) 
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FIGURE 8. If the projectile is a composite particle, it can be excited by the interaction. Therefore, this 
graph will interfere with the double scattering graph in Fig. ^ 



where R A is the nuclear radius. Apparently, this condition is fulfilled for sufficiently 



large values of the projectile momentum p in Eq. (|72[). Thus, as it was first found by 
Gribov in [|TP|], Glauber theory receives important corrections at high energy. 

The condition, Eq. (|73|), which determines, whether Gribov's inelastic shadowing 
becomes relevant, leads us to the important quantum mechanical concept of formation 
time, ox formation length. The formation time is the lifetime of the excitation \n) in Fig. [8| 
in the target rest frame and the formation length is the longitudinal distance over which 
the excited state |n) lives. At high energy, of course, both quantities are identical. The 
formation time/length can be determined in a time-dependent and in a time-independent 
approach. 

In the time dependent formulation, one starts from the energy-time uncertainty rela- 
tion, 

AEAf^l. (74) 
The lifetime of the excitation in rest frame of the projectile is given by 

1 (75) 



f M f -M{ 

In order to obtain the formation time, we have to transform zy to the target rest frame, 
by multiplying xy with the relativistic y-factor, 

t f = At = yr f = tx p (76) 

where 

_ Mr+M- 

M= ; 2 ' ■ (77) 
We finally obtain for the formation time 



(78) 



In the time-independent approach, one starts from the coordinate-momentum uncer- 
tainty relation, 

Ap L Az>l. (79) 

The longitudinal momentum transfer was already given in Eq. [72|. It is calculated in the 
following way, 

/ / M 2 r-Mf 

Ap L = ^Je 2 - Mf - ^E 2 -Mj& J ^ . (80) 

According to the uncertainty relation Eq. ([79]), the excited state lives over the longitudi- 
nal extension 

'f^^MpMf <81) 

As expected, the formation length is identical to the formation time given in Eq. (f78[). 

We see from Eqs. ([78]) and ( |8"T] ) that for large initial projectile momentum, the process 
develops at large longitudinal distances in the target rest frame. At the high center of 
mass energies of RHIC and LHC, the coherence length will be much larger than the 
nuclear radius and all scattering processes will be governed by coherence effects (the 
coherence length becomes as long as several hundreds fm). 



2.2. Elementary hadron-hadron scattering at high energies 

All of the formalism presented above is completely independent of the underlying 
interaction. Before concluding this section, we will briefly discuss the main properties 
of hadron-hadron scattering at high energies. Let us begin by listing some empirical 
facts about hadronic cross sections: 

• Total hadronic cross sections are approximately constant at cm energies of order 
y/s ~ 20GeV and slowly rise, o tot ~ s om , up to the highest energies accessible in 
experiment (Tevatron energy, y/s=l.S TeV). 

• The diffraction cone shrinks as energy increases, indicating that the size of the 
hadron increases with energy. 

• The mean transverse momentum of produced particles is approximately constant 
or increases only slowly with energy, respectively. 

The (approximate) constancy of the total cross section in the framework of QCD 
implies that high energy hadronic scattering is dominated by two gluon exchange []22|], 
see Fig. |9| (left). The two gluon exchange model also yields a purely imaginary forward 
scattering amplitude. In order to explain the increase of the total cross section with 
energy, one has to take the radiation of additional gluons into account, see Fig. [5] (right). 
The probability of gluon emission is proportional to a s y ~ (X s lns, where y is rapidity. 
Thus, each gluon radiation in Fig. |9| (right) contributes a factor Ins to the total cross 
section. Resumming an infinite number of gluon emissions ordered in rapidity, one finds 



FIGURE 9. Double gluon exchange (left), yields an imaginary scattering amplitude and a constant 
cross section. The rise ofhadronic cross sections and the shrinkage of the diffraction cone at high energy 
is due to radiation of additional gluons (right). 



that the total cross section behaves like 
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(82) 



n 



where a<* a s . 

This gluon radiation also explains the shrinkage of the diffraction cone. At high 
energy, the ^-differential cross section in hadronic collisions behaves like 



where 5(5) Ins increases with energy. Such a behavior emerges, if the elastic scattering 
amplitude in impact parameter space is given by 



where the effective hadron radius R\(s) increases as a function of energy. Therefore, the 
shrinkage of the diffraction peak suggests an increase of hadronic sizes with energy. 
In QCD, this can be understood in the following way: Gluons are radiated off the 
projectile with different transverse momenta. As rapidity, or energy, increases, these 
gluons perform a random walk in the impact parameter plane and correspondingly, 
the transverse size of the gluon cloud surrounding the projectile increases. This can 
be regarded as a diffusion process in the impact parameter plane, in which rapidity plays 
the role of time. 

The slow increase of the mean transverse momentum with energy is likely to be 
related to asymptotic freedom. Indeed, at large transverse momentum, p L ^> Aq CD , the 
strong coupling constant becomes small, Ct s {Pi) <C 1, which suppresses the production 
of high p ± particles. 



do 

dt 



do 

dt 



t=0 




(83) 




(84) 



FIGURE 10. At sufficiently high energy, emitted gluons can themselves develop showers. In the squared 
amplitude (right), the gluons combine to ladders, which are denoted by zigzag lines. Fusion of gluon 
ladders is the mechanism behind gluon shadowing. A resummation of fan diagrams like in the left figure, 
corresponds to classical solutions ffify of Reggeon Field Theory /|25|/. 



Eventually, the power-like growth of the total hadronic cross section will violate 
the Froissart-Martin bound P3[], which states that as a consequence of unitarity and 
analyticity, total cross sections cannot rise faster than 

o tot <iCln 2 S , (85) 

where C is a constant. At sufficiently high energy, emitted gluons can develop showers 
themselves, see Fig. [T0|. Due to this process, the projectile sees a reduced gluon density in 
the target and the growth of the cross section is slowed down. This effect, in the squared 
amplitude, realizes a QCD realization of Gribov's inelastic shadowing (see Fig. [IT].) 
Even though such unitarity corrections might already be present in proton-antiproton 
scattering at Tevatron, they will be much more pronounced in nuclear collisions at RHIC. 

As the magnitude of these effects increases with energy and/or atomic number of the 
colliding nuclei, the classification of diagrams in terms of individual nucleon-nucleon 
amplitudes (or parton ladders) rapidly starts to lose sense - the non-linear effects 
become extremely important. The treatment of nuclear interactions in this high-density 
regime will be considered in the following section. 



3. CLASSICAL CHROMODYNAMICS OF RELATIVISTIC 

HEAVY ION COLLISIONS 

3.1. QCD in the classical regime 

Most of the applications of QCD so far have been limited to the short distance 
regime of high momentum transfer, where the theory becomes weakly coupled and 
can be linearized. While this is the only domain where our theoretical tools based on 
perturbation theory are adequate, this is also the domain in which the beautiful non- 
linear structure of QCD does not yet reveal itself fully. On the other hand, as soon as 
we decrease the momentum transfer in a process, the dynamics rapidly becomes non- 
linear, but our understanding is hindered by the large coupling. Being perplexed by this 
problem, one is tempted to dream about an environment in which the coupling is weak, 




FIGURE 11. The place of relativistic heavy ion physics in the study of QCD; the vertical axis is the 
product of atomic numbers of projectile and target, and the horizontal axes are the momentum transfer 
Q 2 and rapidity y — ln(l jx) fx is the Bjorken scaling variable). 



allowing a systematic theoretical treatment, but the fields are strong, revealing the full 
non-linear nature of QCD. We are going to argue now that this environment can be 
created on Earth with the help of relativistic heavy ion colliders. Relativistic heavy ion 
collisions allow to probe QCD in the non-linear regime of high parton density and high 
color field strength, see Fig. [TT|. 

It has been conjectured long time ago that the dynamics of QCD in the high density 
domain may become qualitatively different: in parton language, this is best described in 
terms of parton saturation [ ^7j |28) , [29) ], and in the language of color fields - in terms 
of the classical Chromo-Dynamics [Q; see the lectures [ j3T) ] and [ |32| ] and references 
therein. In this high density regime, the transition amplitudes are dominated not by 
quantum fluctuations, but by the configurations of classical field containing large, ~ 
1/oc,, numbers of gluons. One thus uncovers new non-linear features of QCD, which 
cannot be investigated in the more traditional applications based on the perturbative 
approach. The classical color fields in the initial nuclei (the "color glass condensate" 
[[31]]) can be thought of as either perturbatively generated, or as being a topologically 
non-trivial superposition of the Weizsacker- Williams radiation and the quasi-classical 
vacuum fields [|33L |34], |35|] . 



3.1.1. Geometrical arguments 



Let us consider an external probe J interacting with the nuclear target of atomic 
number A. At small values of Bjorken x, by uncertainty principle the interaction develops 
over large longitudinal distances z ~ l/mx, where m is the nucleon mass. As soon as 
z becomes larger than the nuclear diameter, the probe cannot distinguish between the 
nucleons located on the front and back edges of the nucleus, and all partons within 
the transverse area ~ 1 /Q 2 determined by the momentum transfer Q participate in the 
interaction coherently. The density of partons in the transverse plane is given by 

xGAx,Q 2 ) ,« 

where we have assumed that the nuclear gluon distribution scales with the number of 
nucleons A. The probe interacts with partons with cross section o ~ OC s /Q 2 ', therefore, 
depending on the magnitude of momentum transfer Q, atomic number A, and the value 
of Bjorken x, one may encounter two regimes: 

• op A < 1 - this is a familiar "dilute" regime of incoherent interactions, which is 
well described by the methods of perturbative QCD; 

• op A ^> 1 - in this regime, we deal with a dense parton system. Not only do the 
"leading twist" expressions become inadequate, but also the expansion in higher 
twists, i.e. in multi-parton correlations, breaks down here. 




FIGURE 12. Hard probe interacting with the nuclear target resolves the transverse distance ~ 1 / \[Q 
( Q 2 is the square of the momentum transfer) and, in the target rest frame, the longitudinal distance 
~ 1/ (mx) (m is the nucleon mass and x the Bjorken variable). 



The border between the two regimes can be found from the condition op A ~ 1; it 
determines the critical value of the momentum transfer ("saturation scale"[p7|]) at which 



the partem system becomes to look dense to the probe[]: 

,2 „ xG a( X >0%) 



Q >~ a *^r- (87) 



In this regime, the number of gluons from (|S7|) is given by 

x G A (x,Q 2 )^-^Q 2 s R 2 A , (88) 

where Q 2 R A ~ A. One can see that the number of gluons is proportional to the inverse 
of a s (Q 2 ), and becomes large in the weak coupling regime. In this regime, as we shall 
now discuss, the dynamics is likely to become essentially classical. 



3.1.2. Saturation as the classical limit ofQCD 

Indeed, the condition ( |8~7| ) can be derived in the following, rather general, way. As a 
first step, let us note that the dependence of the action corresponding to the Lagrangian 
(0) on the coupling constant is given by 

4 G£ V G£ V d 4 x. (89) 

r 

Let us now consider a classical configuration of gluon fields; by definition, G^ v in such 
a configuration does not depend on the coupling, and the action is large, S ^>h~. The 
number of quanta in such a configuration is then 



where we re- wrote <E% as a product of four-dimensional action density p 4 and the four- 



dimensional volume V 4 . 

Note that since ( |9l]| ) depends only on the product of the Planck constant/! and the 
coupling g 2 , the classical limit/! — > is indistinguishable from the weak coupling limit 
g 2 — > 0. The weak coupling limit of small g 2 = 4na s therefore corresponds to the semi- 
classical regime. 

The effects of non-linear interactions among the gluons become important when 
dfi^n ~ A 2 (this condition can be made explicitly gauge invariant if we derive it from 

the expansion of a correlation function of gauge-invariant gluon operators, e.g., G 2 ). In 
momentum space, this equality corresponds to 

& 2 ~(V> 2 ~(g 2 ) 1/2 = vp^ (9i) 



1 Note that since ~ A 1 / 3 , which is the length of the target, this expression in the target rest frame can 
also be understood as describing a broadening of the transverse momentum resulting from the multiple 
re-scattering of the probe. 



Q s is the typical value of the gluon momentum below which the interactions become 
essentially non-linear. 

Consider now a nucleus A boosted to a high momentum. By uncertainty principle, the 
gluons with transverse momentum Q s are extended in the longitudinal and proper time 
directions by ~ 1 / Q s ; since the transverse area is 7tR\, the four-volume is V 4 ~ tzR\ / Q^. 
The resulting four-density from (^0|) is then 

»*~ a <fr a <i*t~ Q " (92) 

where at the last stage we have used the non-linearity condition (|9T1), p 4 ~ Q%. It is easy 
to see that ( |9~2| ) coincides with the saturation condition (|87l), since the number of gluons 



in the infinite momentum frame N g ~ xG(x, Q s 

In view of the significance of saturation criterion for the rest of the material in these 
lectures, let us present yet another argument, traditionally followed in the discussion of 
classical limit in electrodynamics P5[]. The energy of the gluon field per unit volume 
is ~ E al . The number of elementary "oscillators of the field", also per unit volume, is 
~ ft) 3 . To get the number of the quanta in the field we have to divide the energy of the 
field by the product of the number of the oscillators ~ ft) 3 and the average energy H ft) of 
the gluon: 

gal 

^~=-t. (93) 
k Hco 4 

The classical approximation holds when N~ ^> 1. Since the energy to of the oscillators 
is related to the time At over which the average energy is computed by ft) ~ I /At, we get 

£ a2 »^4- (94) 

Note that the quantum mechanical uncertainty principle for the energy of the field reads 

E al (o 4 ~H, (95) 

so the condition (Q) indeed defines the quasi-classical limit. 

Since E al is proportional to the action density p 4 , and the typical time is At ~ 
using ( p2| ) we finally get that the classical description applies when 

3.1.3. The absence of mini-jet correlations 

When the occupation numbers of the field become large, the matrix elements of the 
creation and annihilation operators of the gluon field defined by 

Av = Tfa A? +ct AO:*) (97) 
^ ka ka ka ka 

k,a 



become very large, 

"k 

so that one can neglect the unity on the r.h.s. of the commutation relation 



N 7 = (cl c 7 ) > 1, (98) 
ka N ka ka' 



c 7 cl -cl c r =1 (99) 

ka ka ka ka 

and treat these operators as classical c— numbers. 

This observation, often used in condensed matter physics, especially in the theoretical 
treatment of superfluidity, has important consequences for gluon production - in par- 
ticular, it implies that the correlations among the gluons in the saturation region can be 
neglected: 

(A(k l )A(k 2 )...A(k n )) ~ (A(fe 1 ))(A(fe 2 ))...(A(fc n )}. (100) 

Thus, in contrast to the perturbative picture, where the produced mini-jets have strong 
back-to-back correlations, the gluons resulting from the decay of the classical saturated 
field are uncorrelated at k, <^Q S . 



Note that the amplitude with the factorization property fllOOD is called point-like. 
However, the relation ( |100D cannot be exact if we consider the correlations of final- 
state hadrons - the gluon mini-jets cannot transform into hadrons independently. These 
correlations caused by color confinement however affect mainly hadrons with close 
three-momenta, as opposed to the perturbative correlations among mini-jets with the 
opposite three-momenta. 

It will be interesting to explore the consequences of the factorization property of the 
classical gluon field ( |100[ ) for the HBT correlations of final-state hadrons. It is likely 



that the HBT radii in this case reflect the universal color correlations in the hadronization 
process. 

Another interesting property of classical fields follows from the relation 

which determines the fluctuations in the number of produced gluons. We will discuss the 
implications of Eq. fllOip for the multiplicity fluctuations in heavy ion collisions later. 



3.2. Classical QCD in action 

3.2.1. Centrality dependence ofhadron production 

In nuclear collisions, the saturation scale becomes a function of centrality; a generic 
feature of the quasi-classical approach - the proportionality of the number of gluons to 
the inverse of the coupling constant ( P0| ) - thus leads to definite predictions [ ]T7| ] on the 
centrality dependence of multiplicity. 

Let us first present the argument on a qualitative level. At different centralities (de- 
termined by the impact parameter of the collision), the average density of partons (in 
the transverse plane) participating in the collision is very different. This density p is 



proportional to the average length of nuclear material involved in the collision, which in 
turn approximately scales with the power of the number N part of participating nucleons, 

p ~ ^arf ^ ne density of partons defines the value of the saturation scale, and so we 
expect 



-s pan 

The gluon multiplicity is then, as we discussed above, is 

dNg S A Q 2 S 
dr\ ~ a s (Q%Y 



Qs ~ N part- d02) 



(103) 



where S A is the nuclear overlap area, determined by atomic number and the centrality 
of collision. Since S A Q 2 ~ N part by definitions of the transverse density and area, from 
(|HE1) we get 

dN s 

-jr~N part lnN part , (104) 

which shows that the gluon multiplicity shows a logarithmic deviation from the scaling 
in the number of participants. 

To quantify the argument, we need to explicitly evaluate the average density of partons 
at a given centrality. This can be done by using Glauber theory, which allows to evaluate 
the differential cross section of the nucleus-nucleus interactions. The shape of the 
multiplicity distribution at a given (pseudo)rapidity r\ can then be readily obtained by 
using the formulae introduced in section 2: 



da 

dn 



J d 2 b &>(n;b) (l-P Q (b)), (105) 



where P Q (b) * s me probability of no interaction among the nuclei at a given impact 
parameter b: 

P (b) = (l-a NN T AB (b)) AB ; (106) 

G NN is the inelastic nucleon-nucleon cross section, and T AB (b) is the nuclear overlap 
function for the collision of nuclei with atomic numbers A and B; we have used the 
three-parameter Woods-Saxon nuclear density distributions p^]. 
The correlation function £P(n;b) is given by 

^(„;*) = ^=i=exp(-^l£), (107) 
'2%an{b) \ 2an(b) J 



here n{b) is the mean multiplicity at a given impact parameter b; the formulae for 
the number of participants and the number of binary collisions can be found in [Bffl]. 



The parameter a describes the strength of fluctuations; for the classical gluon field, 
as follows from ( |101[ ), a = 1. However, the strength of fluctuations can be changed 



by the subsequent evolution of the system and by hadronization process. Moreover, in 
a real experiment, the strength of fluctuations strongly depends on the acceptance. In 



describing the PHOBOS distribution Q46fl, we have found that the value a = 0.6 fits the 
data well. 



In Fig. [[3], we compare the resulting distributions for two different assumptions about 
the scaling of multiplicity with the number of participants to the PHOBOS experimental 
distribution, measured in the interval 3 < |t]| < 4.5. One can see that almost indepen- 
dently of theoretical assumptions about the dynamics of multiparticle production, the 
data are described quite well. At first this may seem surprising; the reason for this re- 
sult is that at high energies, heavy nuclei are almost completely "black"; unitarity then 
implies that the shape of the cross section is determined almost entirely by the nuclear 
geometry. We can thus use experimental differential cross sections as a reliable handle 
on centrality. This gives us a possibility to compute the dependence of the saturation 
scale on centrality of the collision, and thus to predict the centrality dependence of par- 
ticle multiplicities, shown in Fig. |T4[ (see [ [3"7| ] for details). 
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FIGURE 13. Charged multiplicity distribution at ^/s — 130 A GeV; solid line (histogram) - PHOBOS 
result; dashed line - distribution corresponding to participant scaling (x = 0); dotted line - distribution 
corresponding to the 37% admixture of "hard" component in the multiplicity; see text for details. 




3.2.2. Energy dependence 

Let us now turn to the discussion of energy dependence of hadron production. In 
semi-classical scenario, it is determined by the variation of saturation scale Q s with 
Bjorken x = Q s j \fs. This variation, in turn, is determined by the x— dependence of 
the gluon structure function. In the saturation approach, the gluon distribution is related 
to the saturation scale by Eq.( |87|) . A good description of HERA data is obtained with 
saturation scale Q 2 S = 1 -=-2 GeV^ with W - dependence (W = y/s is the center-of-mass 
energy available in the photon-nucleon system) p8| ] 



Q 2 S « w\ 



(108) 
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FIGURE 14. Centrality dependence of the charged multiplicity per participant pair near TJ = at 
t/s = 130 A GeV; the curves represent the prediction based on the conventional eikonal approach, while 
the diamonds correspond to the high density QCD prediction (see text). The square indicates the pp 
multiplicity. 



where A ~ 0.25 -=-0.3. In spite of significant uncertainties in the determination of the 
gluon structure functions, perhaps even more important is the observation [ |48"1 ] that the 
HERA data exhibit scaling when plotted as a function of variable 



Ql 



(109) 



where the value of A is again within the limits A ~ 0.25 -=-0.3. In high density QCD, this 
scaling is a consequence of the existence of dimensionful scale [ E7L [3H| ] 



Q]{x) = Q 2 (x / x y 



(110) 



Using the value of Qf. ~ 2.05 GeV 2 extracted 0] at y/s = 130 GeV and A = 0.25 @ 
used in [[40"|], equation ( |120[ ) leads to the following approximate formula for the energy 
dependence of charged multiplicity in central Au — Au collisions: 



2 dN, 



ch 



N, 



part 



dr\ 



0.87 



r\<\ 



r s (GeV)\ 
130 J 



0.25 



X 



3.93 + 0.25 In 



(GeV) 



130 



(111) 



At = 130 GeV, we estimate from Eq.flTTTl) 2/N part dN ch /dr\ L, 2 = 3.42 ±0.15, 
to be compared to the average experimental value of 3.37 ±0.12 J4q , ^4] , ^3| , |4^7| ]. At 



y/s = 200 GeV, one gets 3.91 ±0.15, to be compared to the PHOBOS value of 
3.78 ±0.25. Finally, at ^ = 56 GeV, we find 2.62 ±0.15, to be compared to [@§] 
2.47 ±0.25. It is interesting to note that formula fll 1 1| ), when extrapolated to very high 
energies, predicts for the LHC energy a value substantially smaller than found in other 
approaches: 



2 dN 



ch 



N, 



part 



dr\ 



10.8 ±0.5; v^ = 5500 GeV, 



(112) 



T]<1 



corresponding only to a factor of 2.8 increase in multiplicity between the RHIC energy 
of y/s = 200 GeV and the LHC energy of y/s = 5500 GeV (numerical calculations show 
that when normalized to the number of participants, the multiplicity in central Au — Au 
and Pb — Pb systems is almost identical). The energy dependence of charged hadron 
multiplicity per participant pair is shown in Fig.|T5|. 
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FIGURE 15. Energ y dependence of charged multiplicity per participant pair at RHIC energies; solid 
line is the result (111). 



One can also try to extract the value of the exponent A from the energy dependence 
of hadron multiplicity measured by PHOBOS at y/s = 130 GeV and at y/s = 56 GeV; 
this procedure yields A ~ 0.37, which is larger than the value inferred from the HERA 
data (and is very close to the value A ~ 0.38, resulting from the final-state saturation 
calculations [p*D[]). 



3.2.3. Radiating the classical glue 



Let us now proceed to the quantitative calculation of the (pseudo-) rapidity and 
centrality dependences [^]. We need to evaluate the leading tree diagram describing 
emission of gluons on the classical level, see Fig. [H 

A 



4> A (x,."i 2 > 




(^ A (x 2 ,(p t -k t f) 



FIGURE 16. The Mueller diagram for the classical gluon radiation. 

Let us introduce the unintegrated gluon distribution (p A (x,kf) which describes the 
probability to find a gluon with a given x and transverse momentum k t inside the nucleus 
A. As follows from this definition, the unintegrated distribution is related to the gluon 
structure function by 



xG A (x,p?[ 



Pt 



dkf (p A (x,kf 



(113) 



when pj > Qj:, the unintegrated distribution corresponding to the bremsstrahlung radia- 
tion spectrum is 



ot. 1 



(114) 



In the saturation region, the gluon structure function is given by the correspond- 
ing unintegrated gluon distribution has only logarithmic dependence on the transverse 
momentum: 



cp A (x 1 kf)^~-A-kt<Q: 



a. 



t 



(115) 



Note that this "mono-jet" production diagra m ma kes obvious the absence of azimuthal correlations in 
the saturation regime discussed above, see Eq. (100). 



where S A is the nuclear overlap area, determined by the atomic numbers of the colliding 
nuclei and by centrality of the collision. 

The differential cross section of gluon production in a AA collision can now be written 
down as [E7L 0T1] 



E ^ = W^\^ J dk ? a s(PA(xirf)(P A (x2i(p- k )t)i ( 116 ) 

where x l2 = (Pf/\A) ex P( = ' =7 ?)' w i m 1 me (pseudo)rapidity of the produced gluon; the 
running coupling a s has to be evaluated at the scale Q 2 = max{k 2 , (p — k) 2 }. The rapidity 
density is then evaluated from fll 16| ) according to 



" N 1 f B%-\, (117) 



dy qv, I d Pt ( I d 3 p 



AA 

where is the inelastic cross section of nucleus-nucleus interaction. 

Since the rapidity y and Bjorken variable are related by In I fx = y, the x— dependence 
of the gluon structure function translates into the following dependence of the saturation 
scale Q 2 on rapidity: 

Q 2 (s;±y) = Q 2 ( S ;y = 0) exp(±Av). (118) 

As it follows from ( |1 18D , the increase of rapidity at a fixed W = ^fs moves the wave 
function of one of the colliding nuclei deeper into the saturation region, while leading 
to a smaller gluon density in the other, which as a result can be pushed out of the 
saturation domain. Therefore, depending on the value of rapidity, the integration over 
the transverse momentum in Eqs. ( |1 16| ),( pT7| ) can be split in two regions: i) the region 



^qcd < k t < Q smin in which the wave functions are both in the saturation domain; 
and ii) the region A « Q s min < k t < Q s _ max in which the wave function of one of 
the nuclei is in the saturation region and the other one is not. Of course, there is also 
the region of k t > Q Sjmax , which is governed by the usual perturbative dynamics, but 
our assumption here is that the role of these genuine hard processes in the bulk of 
gluon production is relatively small; in the saturation scenario, these processes represent 
quantum fluctuations above the classical background. It is worth commenting that in the 
conventional mini-jet picture, this classical background is absent, and the multi-particle 
production is dominated by perturbative processes. This is the main physical difference 
between the two approaches; for the production of particles with p t » Q s they lead to 
identical results. 

To perform the calculation according to ( |1 17| ),( [lT^ ) away from y = we need also 



to specify the behavior of the gluon structure function at large Bjorken x (and out of 
the saturation region). At x — > 1, this behavior is governed by the QCD counting rules, 
xG(x) ~ (1 — x) 4 , so we adopt the following conventional form: xG(x) ~ x~^ (1 — x) 4 . 
We now have everything at hand to perform the integration over transverse momentum 



in ( |117D , ( |116t ); the result is the following 



dN 2 i / Qsmin 

— = const S A Q s min In — j— 
dy s > mm \A 2 QCD 



X 




Q 



.V.///(tt |y| 



(119) 



= Q 2 s {s;y- 
o 2 



where the constant is energy-independent, S A is the nuclear overlap area, Q s 
0), and Q s min r max ^ are defined as the smaller (larger) values of < ]118| ); at y = 

Qsmax = Q 2 ( s ) = Q 2 ( s o) x ( s / s q)^^ 2 - The first term in the brackets in ( |119D originates 
from the region in which both nuclear wave functions are in the saturation regime; this 
corresponds to the familiar ~ ( 1 / a s ) Q 2 R\ term in the gluon multiplicity. The second 
term comes from the region in which only one of the wave functions is in the saturation 
region. The coefficient 1/2 in front of the second term in square brackets comes from k t 
ordering of gluon momenta in evaluation of the integral of Eq. fll 16| ). 

The formula ( |119D has been derived using the form ( |1 15[ ) for the unintegrated gluon 
distributions. We have checked numerically that the use of more sophisticated functional 
form of (p A taken from the saturation model of Golec-Biernat and Wiisthoff p5| ] in 
Eq.( |l 16Q affects the results only at the level of about 3%. 

Since S A Q 2 ~ N part (recall that Q 2 Aq CD is defined as the density of partons in the 
transverse plane, which is proportional to the density of participants), we can re-write 
( |TT9D in the following final form 



dN 
dy 



cN, 



part 
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Q 



A 2 

QCD , 



X 



X 



l+X\y\ 1 



Q?_ g (l+A/2)[y| 



(120) 



with Q 2 (s) = Q 2 (s Q ) (s/s )^/ 2 . This formula expresses the predictions of high density 
QCD for the energy, centrality, rapidity, and atomic number dependences of hadron 
multiplicities in nuclear collisions in terms of a single scaling function. Once the energy- 
independent constant c ~ 1 and Q 2 (s ) are determined at some energy s Q , Eq. Q120| ) 
contains no free parameters. At y = the expression Q119| ) coincides exactly with the 



one derived in [37], and extends it to describe the rapidity and energy dependences. 



3.2.4. Converting gluons into hadrons 

The distribution ( |120D refers to the radiated gluons, while what is measured in experi- 
ment is, of course, the distribution of final hadrons. We thus have to make an assumption 
about the transformation of gluons into hadrons. The gluon mini-jets are produced with 
a certain virtuality, which changes as the system evolves; the distribution in rapidity is 
thus not preserved. However, in the analysis of jet structure it has been found that the 
angle of the produced gluon is remembered by the resulting hadrons; this property of 
"local parton-hadron duality" (see p3"| ] and references therein) is natural if one assumes 
that the hadronization is a soft process which cannot change the direction of the emitted 



radiation. Instead of the distribution in the angle 6, it is more convenient to use the dis- 
tribution in pseudo-rapidity 77 = — lntan(0 /2). Therefore, before we can compare Q119| ) 
to the data, we have to convert the rapidity distribution ( |120D into the gluon distribution 
in pseudo-rapidity. We will then assume that the gluon and hadron distributions are dual 
to each other in the pseudo-rapidity space. 

To take account of the difference between rapidity y and the measured pseudo-rapidity 
77, we have to multiply fll 19Q by the Jacobian of the y «-> 77 transformation; a simple 
calculation yields 

h{m Pt ;m) = , C ° Sh?? (121) 
^ + sinh 2 77 

where m is the typical mass of the produced particle, and p t is its typical transverse 
momentum. Of course, to plot the distribution Q120P as a function of pseudo-rapidity, 
one also has to express rapidity y in terms of pseudo-rapidity 77 ; this relation is given by 



y(T];p t ;m) = ^ ln 



— 2^ + sinh 77 + sinh 77 



— >p- + sinlr 77 — sinh 77 



(122) 



obviously, h{r\;p t ;m) = dy(r\;p t ;m)/dr\. 

We now have to make an assumption about the typical invariant mass m of the gluon 
mini-jet. Let us estimate it by assuming that the slowest hadron in the mini-jet decay is 
the p-resonance, with energy E p = (m 2 + p pt + p 2 p z ) 1 ^ 2 , where the z axis is pointing 
along the mini-jet momentum. Let us also denote by jc- the fractions of the gluon energy 
q carried by other, fast, i particles in the mini-jet decay. Since the sum of transverse 
(with respect to the mini-jet axis) momenta of mini-jet decay products is equal to zero, 
the mini-jet invariant mass m is given by 



in 



jet 



m 



(Y, x i% + E p) 2 - (L x #z+Pp, 



2j^x.? z -(m P) f-/7p jZ ) =2Q s -m eff , 



(123) 



where m 



p.t 



(m 2 + p 2 p j) l l 2 . In Eq. Q123| ) we used that = 1 and q 



q z = Qs- 



0.5 GeV. 



Taking p p z w p p t w 300 MeV and p mass, we obtain m e ^ 

We thus use the mass m 2 ~ 2Q s m e ^ ^Q s -\ GeV in Eqs.( |121|J122D . Since the typical 
transverse momentum of the produced gluon mini-jet is Q s , we take p t = Q s in ( |121| ). 
The effect of the transformation from rapidity to pseudo-rapidity is the decrease of 
multiplicity at small 77 by about 25 — 30%, leading to the appearance of the w 10% 
dip in the pseudo-rapidity distribution in the vicinity of 77 = 0. We have checked that 
the change in the value of the mini-jet mass by two times affects the Jacobian at central 
pseudo-rapidity to about ~ 10%, leading to ~ 3% effect on the final result. 

The results for the Au —Au collisions at y/s =130 GeV are presented in Figs [FT] and 
T8|. In the calculation, we use the results on the dependence of saturation scale on the 
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FIGURE 17. Centrality dependence of charged hadron production per participant at different pseudo- 
rapidity t] intervals in Au — Au collisions at y/s =130 GeV; from /^9|/, the data are from ffSj]. 



mean number of participants at y/s = 130 GeV from []37|], see Table 2 of that paper. 
The mean number of participants in a given centrality cut is taken from the PHOBOS 
paper p6| ] . One can see that both the centrality dependence and the rapidity dependence 
of the y/s =130 GeV PHOBOS data are well reproduced below r\ ~ ±4. The rapidity 
dependence has been evaluated with A = 0.25, which is within the range A = 0.25 -=-0.3 
inferred from the HERA data [^Sp. The discrepancy above r] ~ ±4 is not surprising 
since our approach does not properly take into account multi-parton correlations which 
are important in the fragmentation region. 

Our predictions for Am — Au collisions at y/s = 200 GeV are presented in [49[]. The 
only parameter which governs the energy dependence is the exponent A, which we 
assume to be A ~ 0.25 as inferred from the HERA data. The absolute prediction for the 
multiplicity, as explained above, bears some uncertainty, but there is a definite feature of 
our scenario which is distinct from other approaches. It is the dependence of multiplicity 
on centrality, which around 77 = is determined solely by the running of the QCD strong 
coupling PTp. As a result, the centrality dependence at y/s = 200 GeV is somewhat less 
steep than at y/s = 130. While the difference in the shape at these two energies is quite 
small, in the perturbative mini-jet picture this slope should increase, reflecting the growth 
of the mini-jet cross section with energy 
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FIGURE 18. Pseudo-rapidity dependence of charged hadron production at different cuts on centrality 
in Au—Au collisions at y/s =130 GeV; from /|£^/, the data are from fit(^j. 



3.2.5. Further tests 

Checking the predictions of the semi-classical approach for the centrality and 
pseudo-rapidity dependence at y/s = 200 GeV is clearly very important. What other 
tests of this picture can one devise? The main feature of the classical emission is that 
it is coherent up to the transverse momenta of about \/2 Q s (about ~ 2 GeV/c for 
central Au — Au collisions). This means that if we look at the centrality dependence of 
particle multiplicities above a certain value of the transverse momentum, say, above 1 
GeV/c, it should be very similar to the dependence without the transverse momentum 
cut-off. On the other hand, in the two-component "soft plus hard" model the cut on 
the transverse momentum would strongly enhance the contribution of hard mini-jet 
production processes, since soft production mechanisms presumably do not contribute 
to particle production at high transverse momenta. Of course, at sufficiently large value 
of the cutoff all of the observed particles will originate from genuine hard processes, and 
the centrality dependence will become steeper, reflecting the scaling with the number 
of collisions. It will be very interesting to explore the transition to this hard scattering 
regime experimentally. 



:dN,/dr| W= 130 GeV 




Another test, already discussed above (see Eq. ( |100| )) is the study of azimuthal corre- 
lations between the produced high p t particles. In the saturation scenario these correla- 
tions should be very small below p t — 2 GeV/c in central collisions. At higher transverse 
momenta, and/or for more peripheral collisions (where the saturation scale is smaller) 
these correlations should be much stronger. 



3.3. Does the vacuum melt? 

The approach described above allows us to estimate the initial energy density of 
partons achieved at RHIC. Indeed, in this approach the formation time of partons is 
T Q ~ l/Qs, and the transverse momenta of partons are about k t ~ Q s . We thus can use 



the Bjorken formula and the set of parameters deduced above to estimate [37] 



for central Au — Au collisions at ^Js =130 GeV. This value is well above the energy 
density needed to induce the QCD phase transition according to the lattice calculations. 
However, the picture of gluon production considered above seems to imply that the 
gluons simply flow from the initial state of the incident nuclei to the final state, where 
they fragment into hadrons, with nothing spectacular happening on the way. In fact, one 
may even wonder if the presence of these gluons modifies at all the structure of the 
physical QCD vacuum. 

To answer this question theoretically, we have to possess some knowledge about the 
non-perturbative vacuum properties. While in general the problem of vacuum structure 
still has not been solved (and this is one of the main reasons for the heavy ion research!), 
we do know one class of vacuum solutions - the instantons. It is thus interesting to 
investigate what happens to the QCD vacuum in the presence of strong external classical 



fields using the example of instantons [J35J]. 



The problem of small instantons in a slowly varying background field was first 
addressed in ^T), |52|] by introducing the effective instanton Lagrangian Uyix) 

L{ ff (x ) = J dpn Q (p)dRexp (^-^ p 2 r[^ v R™' g£ v (* )) (125) 

in which n Q (p) is the instanton size distribution function in the vacuum, Tf^l v is the 

't Hooft symbol in Minkowski space, and R aa ' is the matrix of rotations in color space, 
with dR denoting the averaging over the instanton color orientations. 

The complete field of a single instanton solution could be reconstructed by perturba- 
tively resumming the powers of the effective instanton Lagrangian which corresponds 
to perturbation theory in powers of the instanton size parameter p 2 . In our case here the 
background field arises due to the strong source current 7^. The current can be due to 
a single nucleus, or resulting from the two colliding nuclei. Perturbative resummation 
of powers of the source current term translates itself into resummation of the powers 




FIGURE 19. Distributions of instanton sizes in vacuum for QCD with three light flavors (upper curve) 
versus the distribution of instanton sizes in the saturation environment produced by a collision of two 
identical nuclei for c — 1 (middle curve) and c — 2 In 2 (lower curve) with Q 2 S = 2GeV 2 ; from /pi|/. 



of the classical field parameter oc^A 1 / 3 [T51], |53p. Thus the problem of instantons in the 
background classical gluon field is described by the effective action in Minkowski space 

Sgff = J d 4 x f-^Gj v (*)Gj v (*) + L{ ff (x) + lJ eff (x) + J a ^(x^J . (126) 

The problem thus is clearly formulated; by using an explicit form for the radiated clas- 
sical gluon field, it was possible to demonstrate [ |35| ] that the distribution of instantons 
gets modified from the original vacuum one « (P) to 

^(p)= no(P )exp(- 8g? gg To)2 ), (127) 

where T is the proper time. The result Eq. ( |127| ) shows that large size instantons are 
suppressed by the strong classical fields generated in the nuclear collision (see Fig. |T9|)p|. 
The vacuum does melt! 
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